Dissipative Fluid in Brans Dicke theory and late time acceleration 
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We have investigated the possibility of having a late time accelerated expansion phase for the universe. We 
have used a dissipative fluid in Brans-Dicke(BD) theory for this purpose. The model does not involve any 
potential for the BD scalar field. We have obtained the best fit values for the different parameters in our model 
by comparing our model predictions with SNIa data and the also with the data from the ultra-compact radio 
sources. 
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A number of recent observations [Jj]] suggest that the Q m , the ratio of the matter density (baryonic+dark) to the critical density, 
is significantly less than unity suggesting that either the universe is open or that there is some other sources of this missing energy 
which makes fltotai ~ 1- The recent findings of BOOMERANG experiments [^] strongly suggests the second possibility of a 
flat universe. At the same time, the measurements of the luminosity -redshift relations observed for the 50 newly discovered type 
la supernova with redshift z > 0.35 [Q] predicts that at present the universe is expanding in an accelerated fashion suggesting 

■ the existence of the total negative pressure for the universe. 
pL< ' One of the possibilities is the ACDM model consisting a mixture of vacuum energy or cosmological constant A and cold 

dark matter(CDM). But as the vacuum energy is constant and the matter energy density decreases, it is necessary that their 
ratio must be set to a specific infinitesimally small value (10 -120 ) in the early universe so as to nearly coincide today. This 
is so called "cosmic coincidence" problem. Another possibility is "Quintessence" [0], a dynamical, slowly evolving, spatially 

■ inhomogeneous component of energy density with negative pressure. An example is a time dependent scalar field slowly rolling 
down its potential [|5]. Recently a new form of the quintessence called "tracker field" has been proposed to solve the cosmic 
coincidence problem. It has an equation of motion with an attractor like solution in a sense that for a wide range of initial 
conditions the equation of motion converges to the same solution [|]]. There are a number of quintessence models which have 
been put forward and most of which involve a minimally coupled scalar field with potentials either exponential [Q] or power 
law one [g]. There are also treatments with the nonminimally coupled scalar fields with different type of potentials [|o[. It has 
been shown by Di Pietro and Demaret [^] that for constant scalar field equation of state, which is a good approximation for a 
tracker field solution, the field equations and the conservation equations strongly constrain the scalar potential and the widely 

q-yN used potentials for quintessence such as the inverse power law one, exponential, and the cosine form, are incompatible with 
these constraints. Hence it may be worthwhile to search for a model which will not involve any potential arising from a particle 
£jTj! physics scale. 

Negative pressure can also occur if the CDM fluid is not a perfect fluid but a dissipative one. Recently it has been proposed 
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that the CDM must self interact in order to explain the detailed structure of the galactic halos | |12| |. This self interaction will 
naturally create a viscous pressure whose magnitude will depend on the mean free path of the CDM particles. In a very recent 
work Chimento et.al have shown that a mixture of minimally coupled self interacting scalar field and a perfect fluid is unable to 
drive the accelerated expansion and solve the cosmic coincidence problem at the same time JIT], while a mixture of a dissipative 
CDM with bulk viscosity and a minimally coupled self interacting scalar field can successfully drive an accelerated expansion 
and can solve the cosmic coincidence problem simultaneously. An effective negative pressure in CDM can also be created from 
Cosmic-anti friction which is closely related with the particle production out of gravitational field and one can have a similar 
dynamics like ACDM model as a special case of this Cosmic-antifriction [Q. 

The present work investigates the possibility of obtaining an accelerated universe in Brans-Dicke (BD) theory with a dissi- 
pative fluid. Previously Bartolami et. al. and Sen et. al. [|l0|] have investigated such possibility in BD theory with a perfect 
fluid. But both of them have considered the potential for the BD scalar field itself which was not so in the original BD theory 
JTs]]. But in this work, we have not considered any potential for the BD scalar field. We have compared our solutions with the 
experimental data [^|] to constrain the different parameters in our model. This simple enough model can be useful if one has to 
explain the quintessence model without scalar field potential. 

For a flat FRW universe, with a scale factor R(t), assuming the matter content is a dissipative fluid with only bulk viscosity, 
the BD field equations are: 
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where and p^ are the energy density and pressure corresponding to the BD scalar field and is given by 
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The energy conservation equation for the matter field, which is not an independent equation but can be obtained using (l)-(3) is 
given by 
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We are considering a late time matter dominated universe hence p m = in our case. 

Dissipative effects in FRW cosmology i.e negative tt can be modelled in two ways: Firstly the conventional bulk viscous effect 
in a FRW universe can be modelled within the framework of nonequilibrium thermodynamics proposed by Israel and Stewart 
Jl6|]. In this theory, the transport equation for the bulk viscous pressure tt takes the form 
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where the positive definite quantity £ stands for the coefficient for the bulk viscosity, T is the temperature of the fluid, and r is 
the relaxation time associated with the dissipative effect i.e. the time taken by the system to reach the equilibrium state if the 
dissipative effect is suddenly switched off. Provided the factor in the square bracket is small, one can approximate the equation 
(7) as a simple form 
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which is widely used in the literature. One can also assume that the viscous effects are not so large as observations seems to rule 
out huge entropy productions in large scales [[170. The relation r = can be assumed so as to ensure the viscous signal does 



not exceed the speed of light Jl8[ | and also (tH) 1 = v where v > 1 for a consistent hydrodynamical description for the fluid 
Jl^]. With these assumptions, equation (8) becomes 
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Also as demonstrated in a recent paper by Zimdahl et. al. [[i~3|] one can also 
non conserving interaction inside the matter. This may be due to the particle 
the CDM is not a conventional disspative fluid, but a perfect fluid with varyin 
is an event that occurs in the early universe. But Zimdahl et. al. have shown 
also cause the sufficiently negative tt to violate the strong energy condition. 

In our case, we are not apriori assuming any specific model for this negative 
tt, we have investigated the possibility of having the accelerated phase of the 
the observational estimates. Using (l)-(3) one can write 



6- 



R 



4 = 2.; 



Quo 



R<p 



(9) 



have a negative tt if there exists a particle number 
production out of gravitational field. In this case, 
g particle number. Substantial particle production 
that extremely small particle production rate can 

tt, rather only assuming the existence of a negative 
universe in BD theory, which is comparable with 
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To solve the system of equations we have assumed the following relation between the scale factor R(t) and the BD scalar field 
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where A and a are constants. With (11), equation (10) becomes 
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where [3 = 12 ~ M °; a 6M " . Equation (12) on integration yields 
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where i? and to are positive constants. One can identify to as the present epoch i.e the age of the universe. Now from (11) one 
write 
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where (f>o = ARq. We will assume <fio = 1 without any loss of generality in our subsequent calculations. The solutions for other 
physical quantites now become 
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In these solutions we have essentially three parameters uj, (3 and a, which are related (see just after equation 12). One has to 
ensure that the universe is accelerating i.e < j3 < 1 and also the density parameters for the matter and the scalar field are of 
the same order at present time i.e Qom ~ ^00- These will constrain the different parameters. 

We obtain the best-fit value of (3 by comparing our model predictions with the SNIa data. We use the high-z data of the 
Supernova Cosmology Project (SCP; Perlmutter et al. (1998)) [j^] and the low-z data from Calan-Tololo survey (Hamuy et al. 
1996) for our study. Of the 60 data points, we use a 54 data points for our analysis (Fit C-D of the SCP data; for details of 
the excluded data points see Perlmutter et al. 1998) [j^]. The SNIa data exist up to z ~ 0.9. At larger redshifts, we use the data 
of ultra-compact radio sources (0.55 < z < 3.32; 16 measurements) to constain the the value of (3 [[22|]. 

The x 2 / dof of the comparison of our model with the two data sets is shown in Figure 1. The joint analysis of the two data 
set gives a best fit value (3 — 0.8 with x 2 / dof — 1.18. The good-of-fit probability for the fit Q = 0.12, and the la error on (3 is 
A/3 = 0.05. 

In a very recent paper, using the data for the angular power spectrum of the cosmic microwave background obtained by 
MAXIMA- 1 together with the measurements of high redshift Supernova, Balbi et. al [ pT] | have constrained density parameter for 
matter to be 0.25 < f2 m o < 0.5. We have used this range of f2 OT o together with the value of j3 obtained above by fitting our 
model with different observation, to constrain ui and a. 

We have plotted in figure 2 the parameters u> and a for (3 = 0.5 and (3 = 0.8 and also for Q m o — 0.3 and £l m o = 0.5. For this 
we have used the relation between (3, uj and a and also equation (15). The two values of j3 correspond to the current age of the 
universe 28 Gyr and 18 Gyr respectively where we have assumed that the present Hubble constant is Ho ~ 0.67 x 10~ 10 per year. 
One can see the ranges of the two parameters a and uj, for which these values of (3 and il m o are consistent, are —1.2 < a < —0.8 
and —2.5 < uj < —1.5 approximately. For these ranges of a, the present day variation of G, |^|o = \a\Ho < 10~ 10 per year 
p3|]. One can also check from equations (17) and (18) using these range of parameters that both and 7r remain negative in 
these ranges. Also one can also write from equation (9) 

_ 3[-5 - 5a + a 2 (l + 3w + uj 2 )) 
(uja — 3) (aw + a — 1) 
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FIG. 1. x 2 1 d°f i s shown for SNIa {thin solid line) and Ultra-compact radio sources (dahsed line). The thick solid line shows the result of 
joint analysis of the two data sets 



In Order to have v > 1 which is essential for hydrodynamical description if the CDM is assumed to be a conventional viscous 
fluid, one can not have a particular range for a and u) consistent with the ranges given above. Instead, for a particular value for 
a within the range given above, one can have a range for u>. As an example, for a = —1.2, the range for uj to have v > 1 is 
-2.25 < uj < -1.8. For a = -1, it is -2.5 < u> < -2 and for a = -0.8, it is -2.5 < u> < -2.25. One can see that these 
ranges of uj is consistent with the range shown in figure 2. 
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FIG. 2. The parameter space uj-a for different values of j3 and fl m o\ 1) !2m0 = 0.5, 2) fi m o = 0.3, 3) (3 = 0.8, 4) fj — 0.5 

In conclusion, in recent years it has been shown that a mixture of perfect fluid and quintessence may be an interesting 
candidate to explain a spatially flat universe currently expanding in an accelerated manner. In these models, a minimally coupled 
scalar field rolling down its potential has been used to drive the accelerated expansion and also to account for the missing 
energy of the universe. But all these models necessarily require several fine tuning of different parameters. Also in 
a very recent work, Chimento et.al jTH l have shown that one can not simultaneously solve the cosmic coincidence problem 
and have a late time acceleration in FRW cosmology with a mixture of perfect fluid CDM and a Q-matter. On the other 
hand, recent investigations have predicted that CDM should be self interacting rather than collisionless, in order to successfully 
explain the less dense galactic halos. Hence it is not unreasonable to think that this self interaction may give rise to dissipative 
pressure tt at cosmological scales. In their work, Chimento et.al have shown that a mixture of dissipative CDM and a Q-matter 
can indeed explain the late time acceleration and can solve the cosmic coincidence problem simultaneously. But all of these 
quintessence models whether mixed with perfect fluid or dissipative fluid suffer the problem of unwanted long range forces and 
the quintessence can not be as homogeneous as it should [0]. 

In this work, we have investigated the possibility of having a late time acceleration without any quintessence fields. We have 
used a dissipative CDM model in BD theory for this purpose. The viscous pressure together with negative pressure due to 
the BD scalar field drive the late time acceleration and BD scalar field has been used to account for the missing energy of the 
universe. We also have not used any potential for the BD scalar field unlike the other nonminimally coupled scalar field models 
in literatures JlCj] . The model is simple enough and does not require much fine tuning. We have three arbitary parameters in 
our model which are related through an equation. We have constrained one of the parameters (3 by fitting our model with the 
experimental data from Supernova and also from the ultra compact radio sources. The other two parameters a and uj have been 




constrained using this value of j3 and also assuming that the density parameters due to matter and the BD scalar field are of the 
same order today. These constraints give negative values for us. However, standard limit on us is us > 500 to account the solar 
system tests which sets tight constraints on Post Newtonian deviations from general relativity [^3j : 

1 70 - 1| < 2 x 10~ 3 and \j3 - 1| < 2 x 10~ 3 

where (3$ and 70 denote the usual post Newtonian parameters. But these constraints come from the weak field limit of the theory. 
One should also keep in mind that in extended inflation, the model of La and Steinhardt [ p7[ ] worked provided that us takes a 
value close to 20 which is also not compatible with the solar system tests. It should also be mentioned that in order to explain the 
structure formation successfully in this scalar tensor theory the constraint on us is not at all compatible with the solar system tests 
p^ ] . A negative us is also predicted by the effective models comming from the Kaluza-Klein and superstring theories [ [29| ] . Hence 
it remains always a problem to find a compatibility between the astronomical observations and cosmological requirements. 

The problem is to apply a theory in different scales (Astronomical and cosmological) whereas experiments has been so far 
made only for astronomical scales. We have applied the theory to cosmological scales where still now there is no experimental 
tests for these scalar tensor theories and future datas from supernova at higher redshift may confirm or rule out existence of the 
scalar partner for the graviton. 

It is also important to note that as f3 is constant in our calculation, universe is always in the accelerating phase which seriously 
contradicts with the primeval nucleosynthesis and the structure formation scenario. One way to avoid this problem is to consider 
us not a constant but a function of the scalar field tp. In a recent interesting paper [pOQ, Banerjee and Pavon have shown that with 
us as a polynomial function of cj), one can get both radiation dominated era in the early time and accelerating phase in late time. 
But in that case also asymptotically us acquires a small negative value to have a late time accelerating phase. 

Allowing us to be a function of <f> or redshift z to have both decelerating and accelerating phases at different times while local 
inhomogeneities giving rise to large value of us consistent with solar system test, should be the complete investigation. But this 
will involve a detail computational efforts which is beyond the scope of this paper. What we want to stress that if we can explain 
the quintessence in Brans-Dicke theory without any potential( future observations will predict whether we can or can not) then 
even if us is scale dependent but then also in some scale it has to satisfy the constraints given in our paper in order to explain the 
late time acceleration and also the cosmic coincidence. 
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